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Recent results of N-body simulations have shown that current theoretical models are not able 
to correctly predict the amplitude of the scale-dependent halo bias induced by primordial non- 
Gaussianity, for models going beyond the simplest, local quadratic case. Motivated by these discrep- 
ancies, we carefully examine three theoretical approaches based on (1) the statistics of thresholded re- 
gions, (2) a peak-background split method based on separation of scales, and (3) a peak-background 
split method using the conditional mass function. We first demonstrate that the statistics of thresh- 
^ ■ olded regions, which is shown to be equivalent at leading order to a local bias expansion, cannot 

_ O ' explain the mass-dependent deviation between theory and N-body simulations. In the two formu- 

I lations of the peak-background split on the other hand, we identify an important, but previously 

overlooked, correction to the non-Gaussian bias that strongly depends on halo mass. This new 
term is in general significant for any primordial non-Gaussianity going beyond the simplest local 
/nl model. In a separate paper [1], we compare these new theoretical predictions with N-body 
simulations, showing good agreement for all simulated types of non-Gaussianity. 
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O ' I- INTRODUCTION 

^ , Ongoing and future galaxy surveys will provide a large amount of data that can be exploited to constrain the 

I— I' physics of inflation and the very early Universe, in particular through a measurement of the shape and amplitude 
of primordial non-Gaussianity (NG). Over the past few years, galaxy clustering has emerged as the most powerful 

■ large-scale structure probe of primordial non-Gaussianity (e.g., [2-4]; for a review, see [5]). In particular, references 

] [4, 6-8] have shown that the local quadratic coupling /nl0^ induces a scale-dependent bias 

^: ^^^^ 

o 



in the large-scale power spectrum of biased tracers. Here, bf is the (Eulerian) linear bias factor, Sc ~ 1.69 is the linear 
critical density contrast for spherical collapse, and A4{k,z) cx D{z)k^T{k) is the transfer function between density 
and the gravitational (Bardeen) potential perturbations. Numerical studies have confirmed the scaling Abj cx fc~^ 
and the redshift-dependence Abi cx D{z)~^ [4, 9-13], even though the exact amplitude of the effect remains somewhat 
debatable (at the ^ 10 — 20% level; presumably related to the choice of halo finder [5]). 

However, for other non-Gaussian models such as a local /NLt/*^ model with /c-dependent /nl, the local cubic 
5h ■ coupling g]S!h4>^ , or the orthogonal type, there is a much larger discrepancy between the analytical predictions based 
on the statistics of high-threshold regions [6] and the non-Gaussian bias measured from simulations [14-16]. In the 
5nl</'^ case, the magnitude of the non-Gaussian scale-dependent bias Abi is significantly suppressed relative to the 
theoretical expectation on large scales {k < 0.01 /iMpc~^), even for highly biased halos. The ratio of the measured to 
the predicted non-Gaussian bias strongly depends on the halo mass M: it decreases towards low mass halos, and even 
reverses sign for halos with bf <2 [14]. For the quadratic coupling with fc-dependent /nlI^) oc fc"^', the discrepancy 
between the simulated bias and the high-peak expectation also becomes more severe as the halo mass decreases. 
Furthermore, the deviation depends on the sign and amplitude of the spectral index nf [15]. Recent numerical 
simulations implementing the orthogonal bispectrum shape also show systematic deviations in the measured halo bias 
from the high-peak expectation, in a way that the deviation becomes larger towards lower halo masses [16]. 
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In this paper, we present a careful (re-) derivation of the effect of local and non-local primordial non-Gaussianity 
on the large-scale clustering of tracers (such as galaxies and clusters of galaxies) using the thresholding approach 
[6, 17, 18], as well as two distinct albeit related formulations of the peak-background split (PBS). For all three 
approaches, we present general expressions for the non-Gaussian scale-dependent bias and apply them to models for 
which N-body simulations have been performed. 

In the thresholding approach, we directly calculate the two-point correlation function of halos from the probability 
of finding a single smoothed region above some threshold, and the probability of finding two separate regions above 
the same threshold. This approach has the advantage that the thresholding process is a well-defined mathematical 
operation so that we can in principle calculate the correlation functions without any further approximation. We derive 
a general expression for the amplitude of the non-Gaussian bias in terms of the primordial iV-point functions and the 
Gaussian bias parameters of the thresholded regions, without relying on the high-peak assumption usually assumed in 
previous studies (e.g. [6]). As the deviation between N-body simulations and the theoretical expectation is stronger 
for lower mass halos, such an extension of the high-peak formulation could be seen as a possible resolution. We also 
show that to leading order in /nl, Qnl, ■ ■ ■ , the thresholding approach is equivalent to a local bias expansion. 

In the first PBS approach [4, 7, 19], we decompose the non-Gaussian perturbations into parts that are linear, 
quadratic, and cubic in Gaussian fields. We separate long- from short-wavelength perturbations (the two are uncor- 
related for Gaussian fields, but correlated in the non-Gaussian case), and calculate the bias as the response of the 
halo number density to a long-wavelength density perturbation. This approach is conceptually simple and offers a 
clear physical picture of the impact of primordial non-Gaussianity on the clustering of biased tracers, by isolating 
the effect of the mode-coupling induced by non-Gaussianity. For example, for a generic primordial bispectrum, the 
variance of the small-scale density field is locally rescaled by long-wavelength potential fluctuations. Depending on the 
exact shape of the bispectrum, this rescaling can be scale- independent (local model of non-Gaussianity), which then 
leads to a scale-dependent bias as in Eq. (1); or scale-dependent (e.g. orthogonal and equilateral models), generally 
softening the 1/fc^-dependence. We will see that for cubic-order non-Gaussianity, long wavelength perturbations not 
only rescale the local variance of the density field, but also induce a local skewness. Since the abundance of halos 
also depends on the skewness of the small-scale density field (a fact exploited when searching for non-Gaussianity 
using the mass function of e.g. galaxy clusters), this effect contributes to the non-Gaussian halo bias. This first PBS 
approach has the advantage that it can be generically applied to any prescription for the average halo abundance 
(mass function). On the other hand, it assumes a clear separation between long- and short- wavelength modes, which 
breaks down when measuring the clustering on sufficiently small scales. 

Our second PBS approach is inspired by a calculation of the scale-dependent bias factors in the Gaussian peaks 
model [20] (e.g., the first order bias is bi(fc) = bio + boik'^). In this approach, we apply the peak-background split 
directly to the non-Gaussian density field. This is done by calculating the non-Gaussian conditional mass function 
using an Edgeworth expansion of its Gaussian counterpart. The halo density contrast is then obtained by taking the 
ratio of the unconditional to conditional mass function, and expanding with respect to the large-scale density contrast. 
This allows us to determine the linear bias as the lowest-order coefficient in this series. This approach can in principle 
be applied to any excursion-set mass function. As a first step, we will here formulate it under the assumption that 
the Press-Schechter multiplicity function describes halo abundances. On the other hand, this approach does not rely 
on a separation of scales. Thus, the two PBS approaches presented here make complementary assumptions. 

The paper is organized as follows: we begin by reviewing the models of primordial non-Gaussianity considered here 
and spelling out our notation in Sec. II. We discuss the thresholding approach to non-Gaussian bias and point out 
its limitations in Sec. III. Sec. IV introduces the first PBS approach based on a separation of scales, while Sec. V 
presents the second PBS approach based on conditional mass functions. Sec. VI presents a comparison of the PBS 
and thresholding approaches. The comparison of our predictions with the results of N-body simulations is the subject 
of a companion Letter [1]. We conclude in Sec. VII. 

II. PRELIMINARIES 

A. Four types of primordial non-Gaussianity 

Throughout the paper, we will apply our results to the following four models of primordial non-Gaussianity (NG). 
We parameterize primordial NG via the A^-point functions (A^ > 3) of the Bardeen potential $(a;), a relativistic 
generalization of the Newtonian gravitational potential, in the matter-dominated era. Note that $(a5) has the opposite 
sign relative to the usual Newtonian gravitational potential. Since our goal is to compare analytic predictions with 
the outcome of N-body simulations, our set of models includes all the templates for which simulations have been 
performed. Our main theoretical results, however, will always be given in terms of general A^-point functions and can 
be straightforwardly applied to any given model of non-Gaussianity. 
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1. Local non-Gaussianity 

In local primordial NG, the non-Gaussian field $ is defined by a local Taylor expansion around a Gaussian random 
field <j) as [21-24] 



$(x) = 0(x) + /nl0'(x) +gNL0'(x) . 



(2) 



Here, /nl and (Jnl are dimensionless, phenomenological parameters which we seek to constrain using cosmic microwave 
background (CMB) or large-scale structure (LSS) observations. This type of non-Gaussianity is typically produced 
in inflationary models with more than one scalar field. Since the primeval curvature perturbations are of magnitude 
0(10^^), the cubic order correction is negligibly small compared to the quadratic one when ©(/nl) ^ ©(snl)- 
However, this condition is not satisfied by some multi-field models such as the curvaton scenario, in which a large (/nl 
and a small /nl can be simultaneously produced [25-29]. At leading order, the quadratic term generates a 3-point 
function or bispectrum. 



4=^\ki,k2,k3) = 2/nl P^iki)P4,{k2) + (2 eye.) 



(3) 



where (eye.) denotes cyclic permutations of the indices, Pcf,{k) oc is the power spectrum of the Gaussian 

field 0(x), and rij is its logarithmic slope. On the other hand, the cubic-order terms generate a 4-point function or 
trispectrum. 



^^''\ki,k2,k3,k4) = 65NL P<t>{ki)P^(k2)P^(k3) + (3 eye. 



(4) 



Both bispectrum and trispectrum are peaked on squeezed triangle or quadrilateral configurations, i.e. configurations 
where one side |ki| is much shorter than the other sides. 

2. Scale- dependent /nl 

Next, we will consider a model in which the quadratic coupling dominates but /nl is fc-dependent. The primordial 
bispectrum takes the form [15] 



Cr(kl,k2,k3) = /NL(fcl)P0(fc2)P0(fc3) + (5 perm.) with /NL(fc) = /NL(fcp) 
where kp is some arbitrary fixed scale, and tij is a spectral index. 



(5) 



3. Folded and orthogonal non-Gaussianity 
As a third template, we will consider the folded or fiattened shape, for which the primordial bispectrum reads [30] 

4'\ki,k2,k3) = 6/NL[(P0(fcl)P4fc2) + (2 eye.)) +3(P^(fci)P^(fc2)P4fc3))'^' 

- (P0(fci)i/3p^(fc2)2/3p0(fc3) -f (5 perm.))] . (6) 

The folded shape approximates the non-Gaussianity due to modification of the initial Bunch-Davies vacuum in canon- 
ical single field inflation (the actual 3-point function is not factorizablc) . This template induces a scale-dependent 
bias on large scales with somewhat weaker fc-dcpendence than the local model. [19, 31]. The orthogonal template 
introduced by [32], 

4%i,k2,k3) = 6/NL[-3(P0(fci)P^(fc2) + (eye.)) - 8(P^(fci)P4fc2)P0(fc3))'^' 

+ 3(F4fci)i/3p^(fc2)2/3p^(fc3) + (5 perm.))] , (7) 

gives rise to a similar non-Gaussian halo bias [19], but roughly twice as large in magnitude and opposite in sign (for 
fixed /nl) [69] 



4 



4- Equilateral non-Gaussianity 



Finally, the equilateral type of non-Gaussianity, which arises in inflationary models with higher-derivative operators 
such as the DBl model, is well described by the factorizable form [33] 



(ki,k2,k3) - 6/NL[-(P0(fcl)P0(fc2) + (eye.)) - 2(P4fci)F4fc2)P0(fc3))'^' 

+ {P^{k^f'^P^{k2f'^P^{h) + (5 perm.)) . 



(8) 



It can easily be verified that the signal is largest in the equilateral configurations fci « ^2 ~ fcs, and suppressed in the 
squeezed limit fca ^ fci sa fc2. 



B. From primordial perturbations to galaxies 

In standard CDM cosmologies, galaxies form inside dark matter halos and this introduces a bias between the mass 
and the galaxy distributions [34]. In what follows, we shall adopt a Lagrangian picture. Namely, we express the 
clustering of biased tracers, such as dark matter halos of mass M collapsing at redshift z, in terms of the statistics of 
the initial density perturbation (5ij(k, z) smoothed on a scale R and linearly evolved to redshift z, where R is related 
to M via M = [ATi/ZypR^ . More precisely, 5 is the fractional density perturbation in synchronous gauge. Thus, the 
Poisson equation provides a relationship between (5fl(k, z) and the Bardeen potential ^(x) via 

5ii(k,z) = Al,?(fc,z)$(k) , (9) 

where 

Mn{k,z)=M{k,z)Wn{k) = l^^M^^W^ik) . (10) 

Here, T{k) is the matter transfer function normalized to unity as fc — > 0, g{z) is the linear growth rate of the 
gravitational potential normalized to unity during the matter dominated epoch, and Wn^k) is a (spherically sym- 
metric) window function with characteristic radius R. We will assume a spherical top-hat filter throughout. Note 
also that the matter power spectrum at redshift z is related to the primordial curvature power spectrum through 
P,nik,z)=MHk,z)P^{k). 

Regardless of the initial conditions, we shall denote the Lagrangian bias factors of dark matter halos by bi, bn, 
while Eulerian bias parameters are denoted as bf , etc. Note that these bias parameters are generally scale- 
dependent. The notation bi, 62, bf, etc. will exclusively designate the Gaussian, scale-independent peak-background 
split biases. In the next Section, we will also use the notation ci, C2, . . . for the mass- weighted, cumulative Gaussian 
bias parameters which appear in the thresholding approach. 

We will describe the abundance of halos through their mass function = dn/dAl which we will assume to be of 
the universal form, i.e. 



dlna. 



OM 



d\nM 



(11) 



where /(u) is the multiplicity function and cfom is the RMS density fluctuation on scale M. 

Unless otherwise specified, we shall adopt in all illustrations a flat ACDM cosmology with f7„j = 0.279, h = 0.7, 
and an adiabatic initial perturbations with spectral index Ug ~ 0.96 and amplitude Ag = 7.96 x 10^^° at the pivot 
point fco = 0.02 Mpc^^ (corresponding to a normalization trg ~ 0.81). These values are consistent with the latest 
CMB constraints from WMAP7 [35]. 



III. STATISTICS OF THRESHOLDED REGIONS 



In this Section, we shall present the derivation of the scale-dependent non-Gaussian bias using the statistics of 
regions above threshold [17, 18], without invoking the high threshold (high peak) approximation. Several concepts 
and results introduced in this Section will be employed later in the paper. 
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A. Probability densities 



In the Press-Schechter approach [36], virialized objects are identified with high-density regions in the hnear density 
field. The two-point correlation function of thresholdcd regions, £,y^{r), can be calculated once the probability Pi 
of finding a region whose overdensity is above the threshold 6c ~ 1-69 [37], and the probability P2 of finding two 
such regions separated by a distance r = |x2 — xi|, are known. It is convenient to express the results in terms of the 
significance (peak height) v = Sc/ctq^, where CTos is the r.m.s. variance of the density field smoothed on scale R^- The 
correlation function is then given by [34]: 



P2i> V,^ 

[Pi{> y) 



- 1. 



(12) 



^>i/(r) is commonly interpreted as describing the 2-point correlation of lialos above mass M corresponding to the 
smoothing length i?s. For any non-Gaussian initial density field, P\ and Pi can be expressed in terms of the A^-point 
connected correlation functions as follows [70]: 



Pi(> v) 
P2{> V, r) 



/27r 



dy exp 



E 

.N=Z 



(-1) 



N 



w\ a 



1 

'2^ 



dyi 



dy2 exp 



00 N 

EE 



m dy^ 



w. 



{N,i 



'(r) 



3JV 



For shorthand convenience, we will hereafter omit the explicit z-dependence of (^^(x) (5i^,(x) and Ms{k) 
We have also defined 



(2,m) 



d''™'(r)/ag. (m-i) 

^i'^"^ =0 (to = or 2) , 

.(^■™)=^(^^")(r)/ao^, (iV>2) 



where 



er''")(r) = ( 5,{x,) ■ ■ ■ Ss{x,) S,ix2) ■ ■ ■ S,{X2) ) 



(13) 
(14) 

(15) 
(16) 



m times 



N—m times 



is the A''-point connected correlation function evaluated at two different locations xi and X2. Note that the correlation 
^(JV,o) _ evaluated at zero lag, and that the probability densities Pi and P2 depend explicitly on the 

smoothing scale Rg through the functions wi^'"^^ and the peak height v = Jc/cos- 



B. Bias parameters for a Gaussian density field 

It is instructive to first perform the calculation for Gaussian initial conditions. Later, we shall use the Gaussian 
bias derived in this Section to identify the coefficients of the non-Gaussian scale dependent bias. 



1. Gaussian bias factors from a peak-background split 

When the underlying smoothed density field obeys the Gaussian statistics, the probability Pi of exceeding the 
threshold ly is given by 

Pi(> v) =^= / da;e-^'/2 = - erfc 

In the peak-background split approach, one considers the effect of adding a long- wavelength (background) perturbation 
5i of characteristic wavelength Ri ^ i?s to the small scale density field (peak) 5^ . Assuming that 5i is independent of 
5s- it is clear that adding 5i is equivalent to reducing the threshold v — >■ {5c — 5i)/uqs] thus, Pi(> v,5i), the probability 
Pi in the large-scale overdensity 5i is given by 




Pi{>v,5i) = Pi (>p-^\ 



(18) 
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Wc define the peak-background split bias factors cat as the fractional change of Pi with 5i via 

1 d^Pi{>u,Si) 



CN 



SO that 



CNil^) = - 



Cos 



N 



Pl(> 1^) 

d^[Pi{> ly)] 



Pi{>v) dv^ 



erfc 



'1 e-.V2 



Here, is the Hermitc polynomial defined by 



dx^ 



(19) 



(20) 



(21) 



Note that we adopt the so-called probabilists ' convention for the Hermite polynomials. It is related to the so-called 
physicists ' convention by 



Explicit expressions for the first five Hermitc polynomials are 



Ho{x) 



Hi{x)^x, H2{x)^x^-1, H3{x) 



3x, Hi^x) 



Since v for large we see that in the high-peak limit [v ^ 1), 

CAT w vHn-i{v)I(t^s ~ l^os'-. 



(22) 

- 6a;2 3 . (23) 
(24) 



The one-point probability Eq. (17) and bias parameters Eq. (20) are cumulative. They describe the number density 
and bias of all peaks above the threshold v at fixed smoothing scale Rs- To relate these quantities to the mass 
function and bias of dark matter halos, we follow Press & Schechter [36] and interpret Pi as the fraction of the 
Lagrangian volume occupied by halos of mass exceeding M . Therefore, the halo number density follows upon dividing 
the derivative of Pi w.r.t. mass by M/p, 



dlno-Qs 



d\nM 



(25) 



where the factor of 2 is introduced to account for the fact that regions with 5 < 5c may be embedded in regions with 
5 > 5c on scale > Rs (clouds- in-clouds). Thus, Eq. (25) is of the form Eq. (11) with f{v) = ^2/7ri/exp(— 1/^/2). 
Conversely, integrating Eq. (25) yields 



1 1"^ 

Pi{>iy)^7r- / dM'M'nh{M') . 
Inserting this into Eq. (19), wc find that the c^r are mass-weighted cumulative bias factors. 



CN 



dM'M'nh{M') 



M 



dM'M'nh{M')bN{M') , 



M 



where 



bN{M) 



1 H 



N+l 



(26) 



(27) 



(28) 



are the peak-background split biases derived from the Press-Schechter mass function. Here, vm and ctqa/ denote the 
significance and r.m.s. density fluctuation on the mass scale M. It is only in the high-peak limit (i^ ^ 1) that the 
mass- weighted cumulative bias cn and the bias bN{M) asymptote to the same values [Eq. (24)]. 

So far, we have not yet specified any prescription for how to go from the bias parameter cat to the clustering of 
tracers. This will be elucidated in the next Section, where we calculate the correlation function of thresholded regions 
directly. 
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2. Gaussian bias factors from the correlation of thresholded regions 

111 this Section, we present the calculation of the two-point correlation function £,>u{r) of thresholded regions 
assuming Gaussian initial conditions, and show that the cumulative mass-weighted biases cjy obtained with the peak- 
background split coincide with the bias parameters arising in £,^i,{r). 

Observing that, for Gaussian initial conditions, all the connected correlation functions ^, 
we can express ^>v{r) as 



^^•"^'^ - [P^{> .)? '^TT 



erfc 



dyi 



dV2 exp 



Ur) 



with N >2 vanish, 

p-k{yl+yl) _ 1 . 



Here, £,sij') = Q is the 2-point density correlation smoothed on scale Rs- On employing the definition of Hff{x), 
we can further simplify the double integration as 



dyi 



dy2 exp 



erfc 



+ Ny2N i^w-iMJ e 



(29) 



N=l 



Therefore, we find that the 2-point correlation function of thresholded regions is given by [38] 

1 N 



Mr) = - 



erfc 



An„2N [tlN-i{J^)^ " 



N=l ^^-^Os 



(30) 



Next, on substituting the expression of the cumulative peak-background split bias factors Eq. (20), we can recast the 
peak correlation function into the series 



Mr) = Y.%Ur)V 



N=l 



(31) 



If we compare the expression for ^>u{r) to that obtained from a local bias expansion [39] of the density (5>y of regions 
above threshold, 



s>A-) = Y.%.['^i-)]"^ 



N=l 



(32) 



we see that the coefficient cat is different from the Cat appearing in the correlation function: when calculating ^>,y = 
((5>i,(xi)5,y(> X2)), the coefficient of [Cs('')]^ includes not only c^, but also terms such as CNCN+2mCroT fo'^ positive 
integers m < N/2. This clearly shows that the bias parameters cn from the peak-background split are to be seen 
as "renormalized" bias parameters [40] which take all the higher order moments into account, and thus truly arc the 
coefficients of the observed correlation function of (in this case) thresholded regions. 



C. Two-point correlation function of thresholded regions with non-Gaussianity 

In the presence of primordial non-Gaussianity, all the correlation functions are in principle necessary to 

determine Pi{> v) and P2{> i^,r). Here, we will restrict ourselves to the leading order corrections linear in the 
correlations We derive a general expression for the scale-dependent non-Gaussian bias induced by a primordial 

A^-point function 



1. Relation to local deterministic bias 



First, we show that the leading order contribution to the two point correlation function of thresholded regions, 
which includes terms linear in the connected correlations functions only, is consistent with the result from a 
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local deterministic bias ansatz. Linearizing the exponential factors in Eqs. (13) and (14), we obtain 



(^/2 



1 wl 



E 



P2{> v.r) 



— crfc 
2 



7V=3 
2 



(33) 



V2. 

oo Af-1 



27r E E 



27r 

(7V,m) 



oo (iV.O) 



N=3 



m 



HN-i{v)e 



N=2 m= 



m]{N — to)! 



Hm-i{v)HM-m-i{v)e " , 



(34) 



where we have neglected terms beyond linear order. Thus, the two-point correlation function of thresholded regions 
reads 



n 

oo Af-l 



erfc I —= 



-1-2 oo N-1 



EE 



iN,m) 



N=2 m 



—'^ jn\{N — to) 



mllN — to)! 

W=2m=l ^ ' 

As can easily be seen, a local deterministic mapping 



(35) 



N=0 



(36) 



yields the same result at leading order (the renormalization of the bias parameters cat discussed in the previous Section 
for the Gaussian case will apply at second and higher order). This shows that, at first order in ^, 



the correlation 



of thresholded regions with primordial non-Gaussianity is equivalent to a local deterministic bias relation. Note that, 
for non-Gaussian initial conditions, an effective first-order bias defined through ci oH = \/ C>i^(^)/Cs(^) is generally 
scale-dependent . 



2. Power spectrum of thresholded regions 



We now Fourier-transform Eq. (35), and investigate the separate terms. For simplicity and without loss of generality, 
we will assume that a single non-Gaussian A'^-point function {N > 3) dominates. We then have 

P,Ak) = c^P.(fc)+E Jri^"_";), a"-Hfc) . (37) 

m— 1 ^ ^ 

Here, Ps{k) = (k)P{k) is the matter power spectrum smoothed on scale Rg- Let us consider the term m ~ 1 
first. It is in fact identical to the term to = iV — 1. We have 

?™(fc)= n (y^jCrHk,k:,...,k^_0 (2^)35^(k + ki + ...+k^_i) 

^-2 . 3 V 

= M.,{k)J[ (^y^X,(fc0jx.(9)er(k,ki,...,k^_2,g;X). (38) 

Here, q — — ki — ... — kAr_2 — k, and X is a set of variables characterizing the primordial A^-point function such as 
/nl,5nl, nj depending on the details of the model of non-Gaussianity. In the second line, we have used the fact that 
the matter A^-point function is related to the A^-point function of the potential $ through 

er)(fci,...,fcjv)= |^nA^.(fc.o^ cr^(fci,---,fcjv;x) . (39) 
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Note that the scahng of fs^'^'' in the large-scale limit (fc — > 0) depends on the scaling of ^^"^ in the squeezed limits 
where one argument (fc) is much smaller than the others (fci, . . . , fc7v_2, q)- 

Next, consider the term with m = 2 (again, it is equal to the term m = iV — 2). A similar calculation leads to 



= n (y|^-^.(fcOj-M.(|k-ki|)A^,(g)Cf)(k-ki,ki,...,kAr_2,q;X), (40) 

where now q — — k — k2 — • • • — kAr_2. In the large-scale (small-A:) limit, |k — ki| ^ fc, so that |s^'^'' approaches a 
constant. One can easily verify that this also holds for all m > 3 terms. On large scales, these terms thus all add white- 
noise contributions to the power spectrum of thrcsholdcd regions, and only the terms with m = 1, — 1 contribute to 
the scale-dependent bias. Note, however, that the white-noise corrections which appear for non-Gaussianity of order 
> 4 can be quite significant [14]. 

A general feature of the non-Gaussian power spectrum of tracers in the thresholding approach is now clear: the 
presence of a primordial A^-point function generates a dependence of P>y(fc) on the bias cw-i through l?^'^\k), 
which depends on the scaling of in the squeezed limit. The former has also been pointed out by [41-43], who 
studied the non-Gaussian bias in the local, constant-/NL model. As shown here, this conclusion also holds if we were 
to consider any local bias mapping of the form Eq. (36) (at leading order in the non-Gaussian A^-point function). We 
can then rewrite Eq. (37) as 



P>,(fc) = c?P,(fc) + 2^^i^e™(fc) 



Ps{k) . (41) 



The factor of 2 comes from the sum of the to = 1, A^ — 1 terms, and we have introduced the shape factor 



|n J^MRih)jMR{q)(i''\k,,--- ,kN-2,q,kz;X) , (42) 



where z is some arbitrary unit vector. Noting that P>y = (cf -I- 2ciAci)Ps to leading order in the non-Gaussian 
corrections, we can identify the scale-dependent correction to the linear bias as 



4CA,_1 2 J-i^^(fc) 

(Af-1)!'^"^ Ms{k) 



Aci(fc) - TTr^<^TTrT ■ (43) 



In the rest of this Section, we derive the non-Gaussian correction to the clustering of thresholded regions for the four 
models of primordial NG we consider in this paper. It will prove useful to define general spectral moments through 

<^Ir dkk'^"+'^P^{k)Ml{k) . (44) 

3. Local non-Gaussianity 

For the cubic local model, described by primordial three- and four-point functions [Sec. II A 1], Eq. (41) becomes 



cl+ic^02CTlM:\k)^i^\k,X) + ^alM-\k)^i^Hk,X) 



Ps{k) , (45) 



where -^i'^^(fc, /nl) is precisely equal to the form factor introduced by [6, 31]. Focusing on the quadratic case first, 

(3) 

note that on large scales, J^s — /nl and the power spectrum for ^nl = becomes 

P>Ak) = [cl + 4.hLCiC2alM:\k)] P,{k). (46) 
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From the above equation, it is clear that the scale dependence of the non-Gaussian bias is Ab{k) cx Aij^{k) oc For 
high thresholds ^ 1 in particular, the Gaussian bias parameters cn approach /cr^g so that we can approximate 
the coefficient of the non-Gaussian correction as C2(Jqs — ci5c- Therefore, we recover the expression of [6], 



P>.(fc) "^'bl 



1 + 4/nl ^' 



Ms{k) 



Ps{k) , (47) 



upon replacing ci with hi (i.e., assuming a narrow mass bin). 
For the local ^nl^"^ model, note that 



Ps{k) = Mlik) [1 + 65NLa2] P^{k) (48) 

The matter power spectrum Ps{k) thus contains cr^ = {(fP'), which has a logarithmic divergence for both large and 
small scales [44]. In reality, the finite survey size and the free-streaming scale of dark matter provide low- and high-fc 
cut-offs. In simulations, the finite box size and the resolution provide such cutoffs [14]. 

On large scales, the shape factor J^i'*' generated by the local trispectrum rapidly converges towards 
(3/4).gNLO-§^5'^iloc' where 

sfL^-^ f^MAki)PAki) f^Msik2)P4k2)MA\h + h,\). (49) 

is the skewness parameter of the density field smoothed on scale Rg, (i5s)/(^s)^i ^ local /nl niodel with /nl — 1. 
Therefore, the non-Gaussian contribution to the power spectrum in a pure (7nl model becomes 

APy,{k) ciC3g^i^a^,si%^Ms{k)P^{k) . (50) 

Note that the non-Gaussian bias also has a scale-dependence of k~^. For high peaks ^ 1, ciCaCgs = cfS'^, and 
we recover Eq. (21) of [14] upon replacing ci by bi. In general however, the correct coefficient in the thresholding 
calculation is the third-order bias C3. 



4- Scale- dependent and non-local non-Gaussianity 

For the fc-dcpcndcnt local bispcctrum Eq. (5), the power spectrum of thrcsholdcd regions is 

P>,(fc) = [cl + 4ciC2alM;\k)Ti^\k, hi^,nf)\ P,{k) 

where the redshift independent function J^s (^,/nl5^^/) is computed from Eq. (42) on inserting Eq. (5): 



1 fNhikp) f d^q 



:Ms{q)Ms{\k- q\)P^{q) 



On large scales, the shape factor converges towards 



Tr(3)^i, f N fc->0 1 hhikp) f d^q 2/ NO / N n,- 1 hhikp) 2 
>ik,hL,nf) = -2 -rr^ Ij^^MPMl ' = -2 ^w^^a. 



-'Os kp 



kp 



(51) 



(52) 



(53) 



where (Tqs is the spectral moment evaluated for a = jl. Therefore, the non-Gaussian correction to the peak power 
spectrum becomes 



AP>,(fc) 4ciC2/NL(fcp)fcp "VLA^7'(fc)n(fc) ■ 



(54) 



This result agrees with that of [15, 19] in the high-peak limit only, for which c\C2 ~ (z^/cros)"^- 

Finally, for the folded, orthogonal and equilateral bispectrum shapes, the power spectrum of thrcsholdcd regions is 
also given by Eq. (51), with Tl obtained from an integration over Eq. (6) and Eq. (8), respectively [31]. In the limit 
fc ^ 1, we can set Mgiji) ~ M.g(k^ and, on expanding P^(<i) at second order in fc/fci, we arrive at 



AP>,(fc) = %AciC2hi.<sk~"'^M,{k)P^(k) , 



(55) 



with A=\, a= (us - 4)/6 « -1/2 (folded shape), A = -2, a = {12^ - 4)/6 « -1/2 (orthogonal shape) and A ^ 2, 
a = (tt-s — 4)/3 « —1 (equilateral shape). Again, we recover the high-peak expression [5] if we take the limit v ^ 00. 
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FIG. 1: Ratio of the non-Gaussian correction to the linear bias predicted by the statistics of thresholded regions to that 
obtained in the high-peak hmit. For a non-zero primordial bispectrum (A'' = 3) and trispectrum (A'^ — 4), this ratio is equal to 
b2O'0a/{bi5c) and fe3cros/(fei(5c), respectively. Note that it depends on the order A'^ but not on the specific shape of the primordial 
correlation function. Results are shown at z = as a function of halo mass Af . The Gaussian bias parameters fejv are computed 
from a Sheth-Tormen mass function. 

5. Beyond the high-peak limit 

Fig. 1 shows the ratio of the non-Gaussian linear bias correction arising from the statistics of thresholded regions 
to that obtained in the high-peak approximation. In the literature, cjv is conventionally replaced with bpf, so that 
this ratio becomes ^20'Qg/ (&i^c) for a primordial bispectrum, and b3aQg/{bi6'^) for a primordial trispectrum. Note that 
these quantities do not depend on the shape of the polyspectrum considered. The results are shown at z = using 
the Gaussian bias factors ^at derived from the Sheth-Tormen mass function [45, 46] with p = 0.3 and q ~ 0.75, via 
the PBS approach. These predictions are clearly at odds with the simulation results: firstly, for the local bispectrum 
shape with constant /nl, there is no evidence of a large suppression relative to the high-peak expression [9-11, 47] 
(the correction factor of ~ 0.8 advocated by [10, 11] likely applies for friends-of- friends halo finders solely; see [5] for 
a discussion). Secondly, the simulation studies of [15, 16] unambiguously show that the correction to the high-peak 
expression depends on the shape of the bispectrum. Thirdly, while the suppression seen in Fig. 1 for M > 10^^ Mq//i 
is qualitatively consistent with that measured for the Pnl<^^ model for highly biased halos [14], the sharp upturn below 
^ lO^'^ Mq/H is inconsistent with the findings of [14] at high significance. This appears to exclude the statistics of 
thresholded regions and, more generally, local bias expansions of the form Eq. (36) as a viable framework to calculate 
non-Gaussian bias corrections, at least for realistic halo masses (< 10^^ Mq). We return to these issues in Sec. VI. 

IV. PEAK-BACKGROUND SPLIT: SEPARATION OF SCALES 

In this Section, we present our second derivation of the non-Gaussian, scale-dependent halo bias, based on the peak- 
background split argument. We show that the fact that the cumulants of the density field depend on the smoothing 
scale Rs induces an important and previously overlooked contribution to the non-Gaussian bias correction. 

In this approach, we make a separation of scales and split all perturbations 6, cj), etc. into their long-wavelength 
(subscript "L") and short- wavelength (subscript "S"') pieces, e.g. 

6 = 6l+Ss, (p = (f>L + (l)s, •■• (56) 

Here, short wavelengths signify the scales which impact halo formation (< 10 — 100 Mpc/ft.), while long wavelengths 
correspond to the scales on which we would like to measure the clustering of halos (> 100 Mpc//i). For a Gaussian 
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density field with independent Fourier modes, the L and S pieces are statistically independent. In the presence of non- 
Gaussianity, this is no longer the case. As we will see shortly, it will be convenient to apply the peak-background split 
to the Gaussian primordial curvature perturbation (f). This approach isolates the effect of mode-coupling introduced 
by primordial non-Gaussianity, allowing for direct physical insights. To avoid confusion, we shall denote the physical, 
non-Gaussian density field by (5, to distinguish it from the Gaussian density field 5 related to the Gaussian potential 
</>. 



A. General cubic non-Gaussianity 

Consider the case of weakly non-Gaussian potential perturbations described via non-zero three- and four-point 
functions. We can capture the non-Gaussian corrections by generalizing the cubic local ansatz Eq. (2) in Fourier 
space: 

$(k) = <^(k) + /nL y'|^C.(2)(ki,k2)<^(ki)0(k2)<5z,(k-ki2) 

+ ^^^/^/^ y'^C.(3)(ki,k2,k3)<^(ki)</'(k2)0(k3)<5z?(k-ki23) , (57) 

where ki2... — ki+k2-l-. . . , and the two kernels w^^-* are related to the three- and four-point function, respectively 
[19]. The relation is in general ambiguous, i.e. different kernels can yield the same three- and four-point functions. 
However, the large-scale limit of the non-Gaussian bias depends on the squeezed limit of the iV-point functions, as we 
have seen in Sec. III. In this limit, the kernels w'-^-', w^^-* are unique [71]. 

One possible choice of kernels, which has the nice property (for analytical calculations) of being fully symmetric, is 

(2)n 1 \ 1 ^$ (ki,k2,k3) 
" ^^^'''^^ = 2fe:AP2 + 2perm. ^''^ 

(59) 

65NL P1P2P3 + 3 perm. 

where in the first line, ^3 = |ki2|, while in the second line, ^4 = |ki23|, and Pi = P^{ki). We have pulled out the 
coefficients /nl and gjsjL for convenience. Note that, in general, the four-point function also contains terms of order 
/nL' which we assume to be included in s}^'^ even though we parameterize the amplitude by a single coefficient ^nl- 
Eq. (58) is analogous to the kernel Vl^(ki,k2) defined in [19], and Eq. (59) is the straightforward generalization to 
the cubic case. Note that we define the kernels in terms of (j) here, while they are defined in terms of (/>o(fc) = T{k)(j){k) 
in [19]. The final result (in the large-scale limit) is independent of this choice of kernel, which yields w^^^ = uj^'^^ = 1 
for the local model. 

In the next subsection, we first calculate the effect of long- wavelength perturbations (j)L, on the statistics of 
the small-scale density field 6s- We then derive expressions for the non-Gaussian halo bias for general cubic non- 
Gaussianity. 



1. Effect of long-wavelength perturbations on the density field 

We begin by applying the separation of scales, Eq. (56), to Eq. (57). Clearly, for the quadratic part we will 
obtain the combinations (SS), {SL), and (LL), while the cubic part yields (SSS), (SSL), (SLL), [LLL). The terms 
involving L solely do not influence halo abundance (since they do not contribute significantly to the moments of the 
small-scale density field). The terms involving S'-perturbations only increase the variance, skewness, and kurtosis of 
the small-scale density field. They may thus affect the abundance of halos. However, they do so in a scale-independent 
way and, thereby, induce at most a scale-independent bias correction. Hence, in order to derive the (scale-dependent) 
effect of non-Gaussianity on halo clustering, we only need to retain the mixed terms. 

We now want to derive an expression for the non-Gaussian small-scale density field 5s{k) = M{k)<^s- We obtain it 
by multiplying the short- wavelength part of Eq. (57) by M.{k). Next, we apply a trick, noting that M.{k) oc fc^, and 
P = [ki + (k-ki)]2. Thus, 



M(A:) =M(fci)+Al(]k-ki]) + 0(ki • [k-ki]) 



(60) 
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When inserting Eq. (60) into the first Une of Eq. (57), we see that for the local model, where w*^^^ = 1, the last term in 
Eq. (60) corresponds to the Fourier transform of Vi/)^ • V05 (recall that we are only dealing with mixed terms). When 
averaging over a region where V<f>L is approximately a constant gradient, this term vanishes since (ps is uncorrclated 
with (f)L (see also [8] for a different procedure in the local case). Below we will perform precisely such an averaging 
procedure. A similar reasoning can be applied to the non-local case. Hence we will drop this term and its analogs 
in the cubic part of Eq. (57). Note that we have neglected the fc-dependence of the transfer function here. One can 
circumvent this by defining the kernel in terms of as done in [19]. Eq. (60) and its generalization to several ki 
then lead to 

/j3 j 
_ j^^) [5i(ki)</)5(k - ki) + 5s{^i)M^ - ki)] 

/ffih f rPk 
(2^ y (^^^'nki,k2,k-ki2)[20i5i0s + 0i0i5s + 2</.i<^s^s + 5L</'s'^s] • (61) 

In the second line, we have omitted the arguments of 4>, 5 for brevity (the factors in each product are evaluated at 
ki, k2, and k — ki2, respectively). 

In the presence of non-Gaussianity, the statistical properties of 5s can be derived straightforwardly from Eq. (61) 
by taking advantage of the fact that (ps, are Gaussian fields. We will consider a region of "intermediate" size 
R ^ Rs over which the long- wavelength perturbations can be considered constant. This approximation will break 
down when predicting the clustering on scales which contribute significantly to ctos (see the discussion below and in 
[19]). We then calculate the variance and skewness of Ss in the presence of "external" perturbations (pL, 6l- To 
compute the variance for instance, we calculate ((5s(k)(55(k')) and integrate over k. It is sufficient to consider a single 
(for the quadratic terms) or two independent (for the cubic terms) long-wavelength Fourier modes and, hence, omit 
the integrals over fci, fc2- This is because we will eventually take derivatives with respect to single long- wavelength 
Fourier modes in order to derive the non-Gaussian scale-dependent bias. The variance on scale Rs reads 

= {SsmJsm^r = + 4/nl [M^WLik) + SL{^)al^,{k)] 

+ Qg^i^4>L (ki ) (k2 ) , (ki , k2 ) + 6.9nl [0l (ki )(5l (k2 ) + (ka ) (5l (ki )] (ki , ks ) , (62) 

where indicates an average over a given intermediate region of size R. Note that the terms from quadratic 
non-Gaussianity are linear in 0^,1 while those from cubic non-Gaussianity are quadratic in 0^, 5^. For Eq. (62), 
we have defined the following fc-dependent spectral moments (not to be confounded with Eq. (44), which does not 
depend on k): 

r ^3]„ 

(63) 
(64) 
(65) 
(66) 





1 (27r)3 


J^\k,k,)Ml{k,) P^ih) 




1 (27r)3 


UJ^^Hk,ks)Msiks) P4ks) 


c^L(ki,k2) = 


f d^k. 


oj'^^\-kiM,'i^s)Ml{ks)P4,{ks) 


c^50s(ki,k2) = 


1 (27r)3 


UJ^^\kuk2X.)M,{ks)P4.ih) 



In the following, we will ignore the term 4/NL'5Lcr505 since it only generates a very small (< 10^''/nl) scale- independent 
correction to the halo bias. 

At cubic order in Eq. (61), there are two terms of the type (LSS). These terms indicate that the small-scale 
density acquires a skewness (third moment) which is modulated by long-wavelength perturbations. We will only 
include the effect of the first term, 65nl0l'/'s^5, as the second term proportional to 6l(I>s only produces a scale- 
independent correction to the halo bias. The three-point function of the small-scale density field induced by a single 
long- wavelength perturbation (pLiki) is given by 

((5(k)5(k')<5(k"))fl = 3,9NL<^L(ki) {2TTfSD{k + k' + k") M{k)M{k')M{k") 



X 



where "2 perm" indicates the two cyclic permutations of {k, k,' k"). Recall that the subscript R on the expectation 
value indicates averaging over a region where (pL is approximately constant. In deriving Eq. (67), we have used that 
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k = |k' +k"|, so that A^(fc') + Ai{k") « Ai{k) on large scales. Thus, the three-point function of the small-scale 
density field induced by a long-wavelength perturbation in cubic non-Gaussianity is equivalent to that arising in a 
quadratic model of non-Gaussianity described by the effective three-point function 



4%(k,k',k") = /NL,cff(kO {p''(kz,k',k") +c.(3)(k;,k",k')] P4fc')^0(fc") + 2perm} 



(68) 



where /NL.cff = 3gNL0L(k;). Note that generally depends on the scale ki of the long-wavelength perturbation. 

We can now calculate the skewness parameter of the small-scale non-Gaussian density field, taking out the scaling 
with the long- wavelength mode <pL- 



6 



^^Ms{k^)P4k,) 



(27r3 

XC^(3)(kj^ki,k2)X,(|ki+k2 



lMs{k2)P4,{k2) 



(69) 



(70) 



Here, we have noted that (6^) is already linear in ^nl, so that we can set (J^) — aQ^. 

Summarizing, the effect of long- wavelength modes in general cubic non-Gaussianity is to rescale the local small-scale 
variance of the density field [Eq. (61)], as was discussed for the quadratic case in [4, 7, 19]. This rescaling is linear 
in the long-wavelength modes for the quadratic (/nl) term, and quadratic in 6l, 4'l for the cubic ((/nl) term. The 
terms quadratic in the L-modes induce a non-Gaussian correction to the second order bias bn. We will not consider 
this correction here as it does not significantly impact the halo power spectrum. Furthermore, a long wavelength 
mode in a cubic model also induces a local three-point function (skewness) in the density field [Eq. (67)] : observers in 
a region with (j)^ see a local Universe with an effective quadratic non-Gaussianity described by the "primordial" 
three-point function [Eq. (68)]. 



2. Non-Gaussian corrections to ttie linear bias 



Let us now consider the halo abundance ri/i(x) in some region of size R, with Rg R <^ Ri, and Rg being the 
Lagrangian scale associated with a halo mass M . Throughout, we will assume that hh depends only on the matter 
density pn averaged over R, and the moments of the small-scale fluctuations: ctq^, si'^\ ■ ■ ■ . While, in the Gaussian 
case, a perturbation Sl only changes the average density [pji — ^ P_r(1 + ^l)], it also affects all the cumulants of the 
density field when the initial conditions are non-Gaussian. Applying the chain rule, we find 



bi(fc) ' '^"^ 



fih d5L{k) 



dlnpR dln&os d5L{k) dsi^^ dSL{k) 



(71) 



Here, is the average number density of halos of mass M with non-Gaussian initial conditions, and all derivatives 
are evaluated at Sl = 0. Owing to isotropy, bi(A:) only depends on the magnitude of the fc- vector. The first term in 
Eq. (71) is the usual Gaussian bias bi, while the second and third terms yield the non-Gaussian corrections. Thus, 
the non-Gaussian contribution Abi(A:) to the linear bias bi(A:) arises from the dependence of the halo abundance on 
the variance and skewness of the density field. 

Let us deal with the variance first. As we have seen in the last Section, the change in the variance from cubic 
non-Gaussianity is 0{Sj^). Hence, these terms do not contribute to the linear bias and Eq. (61) gives 



dlncTQs 



dSL{k) 



= 2f^L^^M-\k). (72) 

S, =0 '^Os 



Note that this expression in general depends on the smoothing scale Rg or, equivalently, the halo mass M. 

To proceed further, we will restrict ourselves to the case of a universal mass function for Gaussian initial conditions. 
Therefore, the Gaussian halo number density is given by Eq. (11). Throughout this Section, we will not need to specify 
/(z/) explicitly. The non-Gaussian halo abundance fih will thus depend on &os through the significance = Sc/d'os 
and the Jacobian daos/dlnM. Noting that the Gaussian bias is bi = —6^^dlnf{iy)/d\niy, and taking the derivative of 
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Eq. (72) with respect to InM, we obtain 



dlnaos dSL{k) 



Sl=0 



,(fc) ^ 9lnalsik) _ ^ 



(73) 
(74) 



The second term in the square brackets, 2e„s(fc), has previously been neglected [15, 19]. It vanishes in the scale- 
independent local model, for which w*^^) = 1 and a^is = cos) but is non-zero and generally significant for other 
bispectrum shapes. Physically, this term comes about because a scale-dependent rescaling of the variance [Eq. (72)] 
also changes the significance interval diy that corresponds to a fixed mass interval dM. This in turn affects the 
abundance of halos at a fixed mass and thus contributes to the non-Gaussian bias. The term is absent in the results 
of the thresholding approach (Sec. Ill), since the cumulative two-point correlation C>i,{r) is computed at a fixed 
smoothing scale Rs- We return to this point in Sec. VI. 

- ^ (3) 

In order to derive the effect of cubic non-Gaussianity, we need to determine the dependence of on SI , i.e. the 
effect of a primordial three-point function on the average abundance of halos. Different (albeit related) expressions have 
been proposed for the change in the halo abundance induced by primordial non-Gaussianity [48-55]. For dcfiniteness, 
we will adopt the prescription of [49] derived from an Edgeworth expansion of Pi(> ly) (see also Sec. V), 



V 



fl;,(p,aos,0)i l + iao.S'p) 



[v'' - 'iv 




(75) 



In principle however, any other prescription for the response of halo number counts to a small-scale skewness of the 
density field could be employed here. From Eq. (75), we derive 



6 dlmih 



<^o. dS. 



(3) 



<5i=0 




Sl=0 



1 

V 

V 



(76) 



In the last equality, we have identified the i^-polynomials with the Gaussian peak-background split biases derived 
from the multiplicity function /(z^) = ■\/2/7r v exp(— zy^/2), since our parameterization of Uh in terms of 5*'^' was 
derived within the Press-Schechter formalism [49]. While for high peaks ^ 1, the first term in the last equality will 
dominate, for more abundant halos the second term can contribute significantly. The latter again arises because of 

(3) 

the dependence of Sl^i on the smoothing scale Rg. 
Finally, using Eq. (69) we find 



(3) 



9<^L(fc) 



= 3gNLA^-^(A:)5(3)(fc). 



(77) 



<5i,=0 



Then, using Eq. (71) together with Eqs. (73) and (76), we can assemble the expression for the scale-dependent halo 
bias in a general, cubic order model of non-Gaussianity: 



Ab|"^(fc) =2/NLX-i(fc) 



2£..(fc)] 

■ / 91nSl2(fc) \ , ' 



(78) 



The superscript (k) emphasizes that this correction is fc-dependent, and distinguishes it from a fc-independent non- 
Gaussian bias which we shall denote with a superscript (t). Note that the terms in the first line of Eq. (78) apply for 
any universal mass function prescription. On the other hand, the coefficients in the square brackets of the second line 

will change if a different prescription for dVanhjdSs is adopted. 
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B. Application to models of non-Gaussianity 



Local non-Gaussiamty 



In the local model [Eq. (2)], the kernels Eqs. (58)-(59) are simply tu 



(2) 



(3) 



1. Thus, au 



uos, and the 



skewness s!^J induced by a long- wavelength perturbation becomes S'gYoc' skewness in a local quadratic model 

with /nl = 1 [Eq. (49)]. A more direct way to derive this result is to note that the (LSS) terms in the second line of 
Eq. (61) are obtained by applying the Poisson equation to an effective non-Gaussian potential 



(3) 



</>S 



(t^S + {3gNL4>L)4'l 



(79) 



This relation tells us that, in the presence of cubic local non-Gaussianity, a region with a long-wavelength perturbation 
(j)L looks like a Universe with a local quadratic /nl = 3(7nl0l. 

Since Su^s = 0, the correction to the first order bias Eq. (78) then simplifies to 



Ab|''^(fc) 



1 



M-\k) , 



(80) 



where we have defined 



d\nS 



(3) ^ 
.s,loc 



91ncros 



hi 



(81) 



The term linear in /nl recovers the well-known result for the local quadratic model (this is due to the fact that 
(Tos/cTos = 1 + 2/NLi^i is scale- independent). However, the term linear in t^NL departs from the high-peak expression 
derived in [14] as it includes a correction involving the logarithmic slope of 5*^ on cos- We will return to this point 
in Sec. VI. 



power-law relation, 5*^ 



1 — 10 /i ^Mpc, the scale dependence of S\. {^^ is accurately reproduced by an empirical 
for our fiducial cosmology (this agrees with the findings of [56, 57]). 



In the range Rs 

J^L « 3.08 X 

Hence, the second term in Eq. (81) is approximately 0.145 6i and, therefore, not negligible. 

The ratio of the peak-background split prediction to the high-peak result is given by a-Q^es / {biS"^). In Fig. 2, the 
value of this ratio in the limit fc — >■ is shown as the solid curve. We assume a critical collapse density 5c = 1-69 and, 
in the calculation of the Gaussian biases bN, we employ again a Sheth-Tormen multiplicity function with p = 0.3 and 
q = 0.75. As can be seen, the ratio depends strongly on the halo mass AI. At the redshift assumed here {z = 0), it 
reverses sign around A/ ~ 7 x 10^^ Mq/H. 



2. Scale- dependent and non-local non-Gaussiamty 



We now turn to the other models of primordial non-Gaussianity introduced in Sec. II. Since these are all quadratic 
models, we can ignore the term linear in gNL in Eq. (78). The dependence of Abi(fc) on the shape of non-Gaussianity 
enters through the moment Uujs{k) [Eq. (63)] and through the new correction proportional to d\n(Ji^g{k) / d\n(jQg. Since 
we are interested in the large-scale limit, we can use the scaling of the kernel w*-^' in the squeezed limit in order to 
simplify the analytical expressions. For the local quadratic model with fc-dependent /nl (see Eq. (5)), the kernel in 
this limit reduces to 

a;(k,k,-k)^^°(^^y''. (82) 

As a result, 

<yls{k)='^, (83) 

where a = 7i//2 and the spectral moment (Jas is defined in Eq. (44) with R = Rg. Using the first line of Eq. (78), the 
non-Gaussian, scale-dependent bias correction is then given by 

M-^{k) , a = n//2. (84) 



A6{:l,„Jfc) = 2/NL(fcp)fcp"M — 



bi5c 



2 / dlnags _ ^ 
\ aincTos 
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FIG. 2: Ratio of the non-Gaussian correction to the hnear bias predicted by the peak-background split approach to that 
obtained in the high-peak limit. Results are shown at z = as a function of the halo mass M for a local trispectrum with 
cubic parameter (;nl (solid curve), a local bispectrum with fc-dependent quadratic parameter /nl and index Uf = ±0.6 (dashed 
and dot-dashed curve), the folded and orthogonal template (long-dashed curve) and the equilateral bispectrum shape (dotted 
curve). In contrast to Fig. 1, the ratio sensitively depends on the shape of the primordial A'^-point function. 



The logarithmic derivative of a as w.r.t. cos is always larger (smaller) than unity for a > (a < 0), and reaches unity 
for a 7^ only in the limit ctqs — > 0. For the folded and orthogonal bispectrum shapes [Eq. (6) and Eq. (7)], the kernel 
asymptotes to [19] 



2a 



(85) 



with 2a — (rig — l)/3 — 1 = [ug — 4)/3 w —1, and A — 1 (folded) or A ~ —2 (orthogonal). Note that we have neglected 
corrections of order (k/kg)'^ and higher here (although they are easy to include in a numerical calculation). Inserting 
this result into Eq. (78), we arrive at 



^''5ol/ort(fc)=3A/NL 
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Finally, for the equilateral bispectrum, we have 



(86) 



ci;(k, k^ — k) 3 



with 2a — 2{ns — 4)/3 w —2, which in close analogy with the folded case leads to 

2 



A6i:Ii(fc) = 6/nl ( — 



hi5c 
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a = (n, - 4)/3 . 



(87) 
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Our results agree with those of [15] (for the fc-dependent /nl) and [5, 19] (for the folded and equilateral shapes) apart 
from a factor tal^i^c)^ where 



&i(5c + 2 



\ aincTos 
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(89) 
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This quantity clearly depends on the shape of primordial non-Gaussianity through the slope a. The ratio ea/{biSc), 
which quantifies the deviation from the high-peak approximation, is shown in Fig. 2 for two /e-dependent /nl models 
with Uf = ±0.6 (dashed and dot-dashed curve), for the folded and orthogonal templates (long-dashed curve), and for 
the equilateral (dotted curve) bispectrum shapes. As can be seen, the ratio of the PBS to the high-peak prediction 
depends strongly on bf and the value of a. It is larger (smaller) than unity when a > (a < 0). The suppression 
relative to the high-peak prediction is strongest for the equilateral bispectrum shape, for which a « — 1, but significant 
for all bispectrum shapes we considered. Clearly, this strong mass dependence could be exploited to help constrain 
the shape of the primordial non-Gaussianity. As shown in [1], the results of N-body simulations match the predictions 
derived in this section well. 



V. PEAK-BACKGROUND SPLIT: CONDITIONAL MASS FUNCTION 



In this Section, we consider the third derivation of the non-Gaussian bias based on the conditional halo mass 
function. This is essentially a peak-background split approach since we again consider the effect of adding a background 
perturbation di of characteristic wavelength Ri ^ Rs on the number density fih of biased tracers. However, in 
contrast to the previous formulation, we consider a background density perturbation di which is statistically correlated 
with small-scale density fluctuations. As shown in [20], such a peak-background split approach can be applied to 
obtain the scale-dependent bias factors of (Gaussian) density peaks at all orders. Here, we demonstrate that the the 
implementation of [20] can be generalized to derive the non-Gaussian bias corrections. 

In what follows, fih and Pi will denote non-Gaussian number densities and probability distributions, whereas 
and Pi will designate the Gaussian quantities. Since we will hereafter deal with the non-Gaussian density field, we 
shall revert to the notation of Sec. Ill and simply denote the latter as Ss, Si (and Vs^vi). For simplicity, and since we 
are interested in the behavior on large scales, we shall ignore the peak constraint, which leads to corrections scaling 
as k"^ and higher powers. In other words, we will assume that, for Gaussian initial conditions, the number density of 
virialized objects nhiiy, Rs) identified on the scale Rg follows a Press-Schechtcr mass function. 



A. Cumulants and conditional mass function 



Extending the derivation of the Press-Schechter mass function to the non-Gaussian case, we start from 



hhiiy, Rs) 



p dPi{> iy,M) 



M 



dM 



M dM 



dx Pi{x, Rs 



(90) 



where Pi{v,Rs) is the probability that the linear density contrast of a Lagrangian region of mass M oc R^ equals 
Sc = vctqs, and Pi(> M) is the probability that the same density contrast exceeds Sc. In this Section, we shall use 
the parameter ly = Sc/cros exclusively for the significance corresponding to the critical density with smoothing Rg. On 
the other hand, Vs = Ss/cros and i// = Si/cjqi stand for iV(0, l)-distributed stochastic variables corresponding to density 
perturbations on small and large scales, respectively. As in Sec. Ill, we express the non-Gaussian joint probability 
density P{y) for the A^-dimensional vector of variables y in terms of the corresponding Gaussian probability density, 
by using the following general expansion: 
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where (y^^ • • • y^i,„) ^ are connected cumulants and P{y) is the multivariate Gaussian distribution characterized by the 
covariances (y^iy^^a) [^-g-i 58]. On inserting this expression into Eq. (90), the non-Gaussian mass function becomes 
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In the second line we have assumed that all the cumulants are much smaller than unity. This formula agrees with 
that obtained by [49] at first order. Note that the excursion set approach yields additional, albeit small corrections 
to the Press-Schechter expressions [50]. However, we will ignore them in what follows. 
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We now calculate the conditional mass function nh{i', Rs\i^i, Ri)- By definition, the conditional probabihty for 
having a small scale overdensity Vg on scale Rg given a large-scale overdensity i^i on scale Ri is 

5, Dl n^ P2{Vs,Rs,^l,Rl) ,„„^ 

P{vs,Rs\vi,Ri) = J- — . (93) 

The resulting conditional mass function thus is 

fih{v,Rs\uuRi)^-2j-^— j du,P{us.Rs\yuRi)^-2-^j[PM.Ri)\ J dxP2{x,Rs,i^i,Ri) . (94) 

The joint probability distribution P2{i'^ Rs^i^u Ri) is readily obtained from Eq.(91), 

(00 N / rn N—7n\ \ -f/ \ 



Here, the correlator stands for 

(^s™^^-™), = ^"^o'J"'" (C(^)'5f-'"(a.)), , (96) 

where x is an arbitrary spatial location. The function f{i',i>i,e) is the exponential piece of the Gaussian bivariate 
distribution, whereas Hmn{i^,i^hi) are bivariate Hcrmite polynomials. They can be computed by taking derivatives 
oi f{v,vi,e). Namely, 



We define mixed spectral moments via 



+ vf — 2evui 
2(l-e2) 



1 

al^ = — dk0^+'^P^{k)MMMiik) , 
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quantifying the cross-correlation between small and large scales (the x denotes the splitting of smoothing scales: one 
filter is of size Rs, the other of size Ri). Further, we define the quantity as 

1 

j:1 = — dkk^P^{k)Ms{k)Mi{k)S{k,Rs,Ri) , (99) 
271-^ Jo 

where the form factor S generally is a function of k, Rg and Ri. This definition is broad enough to describe all the 
spectral moments and the cumulants of the density field. For instance, setting S{k, Rs, Ri) ~ k^ yields = ctix ■ In 
the following, we will use the following kernel for E^ : 

S{k,Rs) = A<jl^J'i^\k,X)M-\k) . (100) 
Inserting the expression for the form factor Eq. (42), we see that E^ becomes 

K = <^0s-'<^0l {^^-'yi)c , (101) 

i.e. a mixed A^-th order moment of the density field induced by the primordial iV-point function. 

B. Relative overabundance of rare objects 

The non-Gaussian corrections to the iV-th order Gaussian bias parameters hN can be calculated by expanding 
the relative overabundance of biased tracers nh{i^, Rsli^i, Ri)/nh{t^, Rs) — 1 at order Sj^ . However, throughout the 
remainder of this Section we will consider only the correction to the linear bias. Taking the ratio of the conditional 
mass function to the universal one yields 



^ fihjiy, Rs\vuRi) _ ^ ^ ^ J^dxP2{x,Rs,iyi,Ri) _ ^ 

nh{i^,Rs) PM,Ri)^j;rdxPiix,Rs) 
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Now comes a crucial step in the calculation. As Ri increases, the ratio (v^Vi ~™) /dgj (which is the analog of cr^x /'^o/ 
in the calculation of the peak bias factors) remains finite only if the corresponding form factor S{k, Rs,Ri) does not 
depend on Ri (again, this applies when expanding to linear order in vi). This implies that, in Eq. (102), only the 
terms involving the cumulants {v^) ^ or {y^^^vi) ^ will survive. Therefore, upon taking the limit i?; — > oo, we arrive 
at 



r (X-.. )^ 



V2 



- 1 . (103) 



In order to calculate the non-Gaussian contribution to 6i, it is sufficient to expand the right-hand side of Eq. (103) 
at order bi. The first term appearing in the square brackets can be reexpressed as 
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{v^) HN{x)-NenHN-i{x) + 0{vf) 



In the second line, we successively set e — >■ (we can ignore terms involving e^) and employed the relation H'j^{x) = 
NHn-i(x) to expand the result at first order in vi. To simplify the second term in the curly brackets of Eq. (103), 
we use the fact that /{i', vi, e) in Eq. (97) satisfies the following identity 
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Therefore, 
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We thus obtain 



(iV_ 1)1 \ " '/c " / 

On expanding the numerator of Eq.(103) at first order in vi, we can isolate the Gaussian contribution, which is 
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in agreement with the linear PBS bias for the Press-Schechter mass function derived in Sec. III. We now retain all 
the terms linear in the higher-order cumulants {N > 3) in the linear expansion of Eq. (92) and Eq. (95)) and obtain 
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Using the generating function exp(xf:— <2/2) = HN{x)t^ /Nl, we can easily evaluate the integrals over the Hermite 
polynomials. In particular, we find for N > 2: 
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On inserting this expression into Eq.(109), taking the derivative with respect to M and employing the recurrence 
relation Hfqj^i[x) = xHj^[x) — NHj^_i{x), the conditional overabundance of halos simplifies to 



5h{5i) 



biSi 



^1 E ^hl^c-M-f)>A^-lM + (-i')e^A^M 



N^3 



oo ^ 



V 



(111) 



where a primed variable X' now designates dXjda^s (we have used the fact that dv ~ — Sedans/ (^os 



C. Non-Gaussian bias corrections 



In order to calculate the non-Gaussian bias corrections, we have to compute the derivative of the A^-point cumulants 
{v^)^ ^ ag-^si^'' and (lyf-^iyi)^ = / (cr^-Vo;) with respect to ctq.- These are 



{N-2) + 



(N) 



I N-l \' 1 
{''s ''I), = -N. 



N 



(112) 



Replacing the Hermite polynomials with the Gaussian peak-background split biases inferred from the Press-Schechter 
multiplicity function [Eq. (28)], 
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the conditional overabundance of halos can be recast into 
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We can now read off the scale- independent correction Abj'^ (involving the terms proportional to CTq^/ctq;) and a 

scale-dependent correction Ah\^^ (involving the terms E^/ctq; and (E^)'/(Tqj) to the first order Gaussian bias &i. 

The non-Gaussian bias contribution thus is Abj = Ah^^^ + Abj'^''. After some manipulation, the scale-independent 
non-Gaussian bias correction reads 
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where X is again a vector of variables describing the amplitude and shape of the primordial A^-point function. In 
order to write down an explicit expression for the scale-dependent, non-Gaussian bias correction, we use the definition 
of the kernel S{k,Rs) [Eq. (100)], yielding 
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k [h/Mpc] 

FIG. 3; A comparison between the non-Gaussian scale-dependent bias correction Eq. (118) and its low-A: limit Eq. (78) for 
some of the bispectrum shapes and the local trispectrum considered in this work. In all cases, a dotted curve represents the 
low-fc limit. Results are shown as a function of k for halos of mass M = 5.3 X 10^^ Mq/H at z = 0.5, assuming /nl = 100 (for 
the bispectra) and gNL = 10® (for the local trispectrum). 



Then, by definition of the linear halo bias, (5ft,(k) = b{k)5s{'k), correlating the last line of Eq. (114) with 6i yields 
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from which we can read off the scale-dependent non-Gaussian bias correction as 
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This is the main result of this Section. In the high-peak limit, 67V-2 S> &Ar-3 and the first term in the curly bracket 
dominates. Therefore, we exactly recover the results of [6, 14, 15, 31] for the constant /nl, constant ^nl, folded shape 
and fc-dependent /nl, respectively. The second term in the curly brackets arises owing to the mass-dependence of the 
reduced cumulants si^\ As we will see shortly, this term agrees with the correction derived in Sec. IV in the limit 
fc -> 0. 

Note that [59] also employed the bivariate Edgeworth expansion to explore the effect of a local primordial trispectrum 
on the (configuration space) bias of tracers. However, they did not derive any explicit expression for the non-Gaussian 
bias. 



D. Comparison of the PBS approaches 



Interestingly, if we ignore the mass-dependence of the cumulants, then the /c-dependence of Eq. (118) is exactly 
the same as that predicted by the correlation of thresholded regions (see Sec. III). This follows from expanding the 
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non-Gaussian density field in cumulants, which is also done in the expansion of the correlation function of thresholded 
regions. By contrast, our first formulation of the peak-background split (see Sec. IV) leads to a different fc-dependence 
on smaller scales. This difference arises because we have assumed that the long-wavelength perturbation is constant 
over some intermediate scale R over which the halo abundance is averaged. This is a sensible assumption as long as 
the scale "L" over which we measure the clustering of halos is much larger than the scales that contribute to ctos- 
Then, the kernel lj in Eq. (63) is indeed evaluated in the squeezed limit, ks 3> k, and both PBS formulations agree 
exactly. To see this explicitly, we write Eq. (118) for the cases of iV = 3 and = 4: 

Ah\^\k,Rs,N^-3)^2Ls, + ^^^^\Ti'\k)M-'ik) (119) 



Ah[^\k,Rs,N^A)^^h2S, + bi 
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In the large-scale limit, we can use the same approximations made in Sec. IV, i.e. assume that k is much smaller 
than the scales which contribute significantly to the integrand in Eq. (42). On inserting the definition of the kernels 
a;(2), a;(3) [Eqs. (58)-(59)], we obtain 

-Ff\fc)'-5"/NL^^ (121) 
mk)'=\9^i.olS^^}{k). (122) 

Substituting these expressions into Eqs. (119)-(120), we eventually recover Eq. (78) in Sec. IV. 

On smaller scales k > 0.02/i/Mpc around which the matter power spectrum peaks, the separation of scales "L" 
and "S"' is no longer accurate and the predictions of Eq. (118) diverge from the fc — s- limit. In Fig. 3, the exact 
fc-dependence of the non-Gaussian bias correction predicted by the correlated PBS approach, Eq. (118), is compared 
to that predicted by the low-fc expression, Eq. (78). We can see that the latter is accurate to a few percent at 
wavenumber fc < 0.01 hMpc~^. Only for the folded and equilateral shape does the low-fc expression yield a noticeably 
larger non-Gaussian bias correction on scales fc > 0.01 hMpc^^ . This is also true for the orthogonal template (not 
shown in the figure since it is essentially equal to the folded case). The exact difference, however, depends somewhat 
on halo mass and redshift. A quantitative comparison of the scale-dependent bias predicted by the uncorrelated PBS 
approach with that obtained from the statistics of thresholded regions can also be found in [19] (note however that 
the new term derived in this work is not included there). 

Finally, while in the limit fc — >■ Eq. (118) reproduces the well-known result for the local scale- independent 
/nl nrodel [4, 6], at finite fc this expression receives a negative correction from the second term proportional to 
SlnJ^i'^'' /91ncros that increases with wavenumber. At fc = 0.05 hMj>c~^ for instance, the suppression is ^1% and ^4% 
for biased tracers with bf ^ 2 and 3.5, respectively. 



VI. PEAK-BACKGROUND SPLIT VS THRESHOLDING 

We now compare our final result Eq. (118), with the result from thresholding in the high-peak limit, 

Ab^'^PVfc) - -^^^^a^ -^^"^^W (123) 
(A^-l)!''°^ Ms{k) ' ^ ^' 

obtained from Eq. (43) by replacing cat with 6jv. We see two important differences. Firstly, in the thresholding 
approach (which is equivalent to local biasing), the correction to the halo power spectrum induced by a primordial 
A^-point function is proportional to Bn-i- In the PBS approach on the other hand, the correction comes in through 
the dependence of the halo mass function on the {N — l)-th moment si^~^^ of the small-scale density field. The 
latter is proportional to b]\[-2 when the Edgeworth approximation method is applied to the halo mass function. The 
simulation results for all types of primordial non-Gaussianity simulated so far clearly follow the dependence on 
rather than b^-i, thus favoring the interpretation provided by the PBS approach. 

Secondly, the term proportional to dlnJ^s /dlnaos in the PBS prediction [Eq. (118)] is absent in the thresholding 
approach. In [1], we show that the inclusion of this term yields a good match to the simulated halo bias in non-Gaussian 
models beyond the simplest, local quadratic non-Gaussianity with scale- independent /nl- In the thresholding approach 
on the other hand, we associate the correlation of regions above a threshold Sc{z) in the linear density field smoothed 
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at a fixed scale Rg with that of halos above a mass threshold M{Rs) at redshift z. However, halos spanning some 
mass interval should be identified with Lagrangian regions spanning a range of smoothing scales. Consequently, the 
abundance of halos in a mass interval [M,M + dM] not only depends on the cumulants si^ of the density field 
smoothed on scale Rs , but also on the variation of these cumulants with Rs (parameterized through dlnSs^ /dlnaos ) ■ 
An alternative way of seeing this is to describe the abundance of halos in a non-Gaussian density field through an 
effective significance {'{v, Ss"^^), which is defined upon requiring 

- P dlni/ p dhiC> 

n, = M2/ng(-)^ = Jpfoi-)^ ■ (124) 
In the case of a Press-Schechter mass function [Eq. (92)], iy is given by 

---i^-E^Si-^-l^^^'b^-^] (125) 
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The Jacobian (ilni>/(ilnAf in Eq. (124) involves dlnSi"^^ / dlnv , showing that Uh depends on the scale-dependence of 
the cumulants. 



VII. CONCLUSION 



We have carefully re-examined the derivation of the effect of primordial non-Gaussianity on the large-scale clustering 
of tracers (such as galaxies and clusters) beyond the local /nl model, using the statistics of thresholded regions as well 
as two formulations of the peak-background split (PBS). We have shown that the thresholding approach is equivalent 
to local biasing, when considering the leading order contributions from non-Gaussianity. This approach predicts the 
same scale-dependence as the PBS approach in the limit fc — ^ 0, in agreement with the findings of [19]. However, 
unless we consider the high-peak limit (which is never attained by real tracers), PBS and thresholding predict different 
amplitudes for the non- Gaussian contribution to the linear bias. 

We have presented two complementary peak-background split derivations of the effect of non-Gaussianity. In the 
first approach, the separation of scales is invoked to split the Gaussian density field into uncorrelated short- and 
long-wavelength perturbations. This allows us to isolate the mode-coupling effect responsible for the scale-dependent 
bias induced by non-Gaussianity. In the second approach, the separation of scales is invoked to expand the ratio of 
the unconditional to conditional mass function in terms of large-scale perturbations in the non-Gaussian density field. 
Notice that no assumption of separation of scales is made in the thresholding approach, where biasing is a function of 
the local density only. While in the second PBS approach we have restricted ourselves to the case of a Press-Schechter 
mass function, we have nonetheless been able to identify the non-Gaussian bias correction to the linear bias. Both 
PBS approaches predict exactly the same correction in the limit fc — > (once the Press-Schechter expressions for the 
Gaussian biases are identified with fe^r). While they depart at higher wavenumbers (fc > 0.02 h/Mpc), this deviation is 
not very significant for the local or folded type of non-Gaussianity where the non-Gaussian bias correction is strongly 
suppressed at small scales. 

In both approaches, we uncover a new term depending on the scale-dependence of the small-scale moments of the 
density field induced by non-Gaussianity. Physically, this term is induced by the mapping from local significance 
i> — dc/cTQs to mass M: a scale-dependent modulation of ctos changes the interval dv corresponding to a fixed mass 
interval dM. This correction to the high-peak expression of the linear non-Gaussian bias has not been pointed 
out in any previous work. It can be very large for all the models considered here, except for the local bispectrum 
with constant (i.e., fc-independent) /nl- Moreover, we have found very good overall agreement between the PBS 
predictions and the simulated non-Gaussian halo bias [14-16] for the local gj^^cj)^ model, the local ff^hcj)^ model with 
fc-dependent /nl, and the orthogonal bispectrum. This comparison is detailed in a companion Letter to this paper 
[1]. Consequently, the simulation results rule out thresholding, and more generally local biasing, as a viable approach 
to predicting the impact of primordial non-Gaussianity on halo clustering. These new accurate predictions can be 
combined with optimal weighting schemes [60-65] in order to extract information on the scale-dependent bias from 
numerical simulations and forthcoming galaxy surveys. 

In order to further test the PBS approach with numerical simulations, it will be important to take into account 
the scale-independent correction Abj'^ induced by non-Gaussianity through its impact on the abundance of halos. In 
the case of local cubic non-Gaussianity, it will also be necessary to measure the Gaussian second-order bias factor 62 
directly from the simulations, as the effect on the linear bias scales with 62- 

Finally, a natural generalization of the conditional mass function approach discussed in Sec. V is a derivation of the 
non-Gaussian bias factors within the excursion set formalism, for generic moving barriers and non-Gaussian initial 
conditions [50, 51, 66, 67]. We leave these issues for more detailed future treatments. 
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